混凝土破裂过程的三维数值模型

利用混合同余法产生了随机变量,将混凝土中的骨料简化为球体,采用循环比较法将随机变量赋给球心坐标(X,Y,Z),实现了混凝土骨料的空间随机分布,最终建立了三维混凝土数值模型。对模型整体进行单元剖分,将骨料、砂浆、界面层分别投影到该有限元网格中,采用多次网格划分技术完成模型网格划分。通过对算例进行有限元计算,结果表明该模型基本上反映了混凝土骨料分布的实际情况,较好地模拟了混凝土的不均匀性与各向异性,验证了建立该模型的方法是可行的。     

隐式Newmark法分析波动问题精度的探讨

本文研究了用平均加速度的Newmark积分格式求解一维平面弹性波动问题的精度。由于平均加速度的Newmark积分格式没有振幅衰减、不产生数值阻尼,因此采用此法分析波动问题会产生较强的高频振荡,尤其是加速度的计算值误差甚大。本文对防护工程中常见的平台型荷载和升压三角形荷载作用下的波动问题,应用有限元分析的实例表明:采用有数值阻尼的Newmark积分格式是改进有限元求解波动问题的一个有效方法,可以大大地提高应力和加速度求解精度。文中推荐的有数值阻尼Newmark积分格式的积分参数(=0.9,=0.49)可供工程中土与结构动力分析时参考使用。     

Averaging methods for numerical simulations of flows through heterogeneous porous media

Many natural rock systems contain small patches of different permeability which affect the flow of fluids through them. As these heterogeneities become smaller and more numerous, they become harder to model numerically. We consider how to reduce the computational effort required in simulations by incorporating their effects in the boundary conditions at the edges of each grid block. This is in contrast with current methods which involve often arbitrary changes in the fluid properties. The method is restricted to the case of widely-spaced patches, which simplifies interaction effects. The system then reduces to an array of dipoles, and two averaging methods are proposed for finite grid blocks. Several infinite systems, including vertical and horizontal bands, are also considered as further approximations. There is a great wealth of existing results from different fields which lead to identical mathematical problems and which can be used in these cases. Finally, we consider how to use these techniques when the precise configuration of the grid block is not known, but only its statistical properties. This can lead to results which are very different from the deterministic case.     

The low Reynolds number turbulent flow and mixing in a confined impinging jet reactor

Turbulent mixing takes an important role in chemical engineering, especially when the chemical reaction is fast compared to the mixing time. In this context a detailed knowledge of the flow field, the distribution of turbulent kinetic energy (TKE) and its dissipation rate is important, as these quantities are used for many mixing models. For this reason we conduct a direct numerical simulation (DNS) of a confined impinging jet reactor (CIJR) at Re = 500 and Sc = 1. The data is compared with particle image velocimetry (PIV) measurements and the basic flow features match between simulation and experiment. The DNS data is analysed and it is shown that the flow is dominated by a stable vortex in the main mixing duct. High intensities of turbulent kinetic energy and dissipation are found in the impingement zone which decrease rapidly towards the exit of the CIJR. In the whole CIJR the turbulence is not in equilibrium. The strong mixing in the impingement zone leads to a rapid development of a monomodal PDF. Due to the special properties of the flow field, a bimodal PDF is generated in cross-sections downstream the impingement zone, that slowly relaxes under relaminarising conditions. The time required for meso-mixing is dominating the overall mixing performance.     

生物瓣膜应力的数值模拟

应用有限元方法分析了瓣架的材料特性、瓣叶的材料特性和瓣叶的钙化等因素对瓣叶的应力分布的影响。分析结果表明：采用全弹性材料的瓣架将改善瓣叶的应力分布，最大应力区域将远离接缝部位；将生物瓣瓣叶材料简化为拟线弹性材料时，对瓣叶最大主应力等值线的分布影响不大，计算结果在一定范围内存在差异；当瓣叶出现钙化时，钙化点及其周围产生应力集中，促进钙化的进一步加剧。     

爆轰波拐角传播三维数值模拟

为了研究爆轰波在拐角中的传播特性以及拐角爆轰低压流场现象,用LSDYNA 3D程序对三种常见的拐角装药的爆轰波传播特性进行了数值模拟,清晰显示了爆轰波因拐角处爆轰波传播面积的变化而产生的衰减 增长过程,讨论了侧向三通中爆速的变化情况,并得到了拐角处低压爆轰区的尺寸。     

NUMERICAL SIMULATION OF COMPLICATED PIPE FLOW WITH k-εMODEL

In this paper, turbulence in a complicated pipe is simulated by using the k-ε model. The ladder-like mesh approximation is used to solve the problem of complicated boundary with the result of numerical simulation favorable. Two computational examples are given to validate the strong adaptability and stability of k-ε model.     

超级有限元法及其在结构工程中的应用

本文探讨一种基于半连续半离散思想，适用于复杂结构（如高层框架、剪力墙、桁架、网架等结构系统）工程分析的超级有限元，其结构数值分析是按连续体进行，但又按单个构件进行有限元计算。这种按整体系统进行离散所获得的单元内部包含众多构件，有别于一般常见的实体有限元，称为“超级有限元”。这种方法自由度数比一般有限元法少很多，又与单元内部所含构件数多少无关，并可求取结构内每个构件的内力值。     

推进－迭代算法计算效率的研究

本文对求解三维定常超音速动性流场的一次空间推进，在每一个推进站沿伪时间层局部迭代的推进－迭代算法作了进一步的研究．在每一推进站（侧向平面）沿伪时间层局部迭代时，给出了四种不同的隐式迭代方法，即沿侧面两个方向（法向和周向）全用隐式；法向隐式而周向采用Ｇａｕｓｓ－Ｓｉｌｄｌｅ来回扫描迭代；法向隐式而周向显式及以系数矩阵谱半径代替系数矩阵的简化标量隐式算法．用这四种算法模拟了三维球锥黏性绕流，给出了四种不同算法的计算效率和收敛特性比较．     

FINITE ELEMENTAN ALYSIS FOR THE UNSTEADY NEARSHORECIR CULATION DUE TO WAVE－CURRENT INTER－ACTION（I）──NUMERICAL MODEL

ＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＦＯＲＴＨＥＵＮＳＴＥＡＤＹＮＥＡＲＳＨＯＲＥＣＩＲＣＵＬＡＴＩＯＮＤＵＥＴＯＷＡＶＥ－ＣＵＲＲＥＮＴＩＮＴＥＲ－ＡＣＴＩＯＮ（Ｉ）——ＮＵＭＥＲＩＣＡＬＭＯＤＥＬＷｕＷｅｉ－ｘｉｏｎｇ（吴伟雄）（Ｔｏｎｇｊ...